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Abstract 

Recent combined analyses of the CMB and galaxy cluster data reveal unexpectedly large and 
anisotropic peculiar velocity fields at large scales. We study cosmic models with included vorticity, 
acceleration and total angular momentum of the Universe in order to understand the phenomenon. 
The Zeljdovic model is used to mimic the low redshift evolution of the angular momentum. Solving 
coupled evolution equations of the second kind for density-contrast in corrected Ellis-Bruni covari- 
ant and gauge-invariant formalism one can properly normalize and evaluate integrated Sachs- Wolfe 
effect and peculiar velocity field. The theoretical results compared to the observations favour a 
much larger matter content of the Universe than that of the concordance model. Large-scale flows 
appear anisotropic with dominant components placed in the plane perpendicular to the axis of 
vorticity (rotation). The integrated Sachs- Wolfe term has negative contribution to the CMB fluc- 
tuations for the negative cosmological constant and it can explain the observed small power of the 
CMB TT spectrum at large scales. The rate of the expansion of the Universe can be substantially 
affected by the angular momentum if its magnitude is large enough. 
PACS numbers: 98.80.Es: 12.10.Dm: 04.60.-m 
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I. INTRODUCTION AND MOTIVATION 



The phenomenology of particle physics, cosmology and astroparticle physics should be 
understood with the unique physical principles. Heavy Majorana neutrinos within the non- 
singular and causal SU(3) electroweak-strong unification []| represent the main building 
block of the Universe as a cold dark matter (DM) particles. They are cosmologically stable 
T/Vi >> Tu with large annihilation cross sections [2]. The H.E.S.S. source J1745-290 sl in 
the center of our galaxy, WMAP haze [J] and EGRET excess of diffuse photon background 
are possible indications of the heavy DM particle annihilations. 

Light Majorana neutrinos trigger vorticity of the Universe 5|, IgJ through the spin-torsion 
connection in the Einstein-Cartan cosmology at early times of the evolution. The formation 
of large-scale structures in the form of galaxies and clusters, as well as, the present anisotropy 
of spacetime, indicate the existence of nonvanishing total angular momentum of the Universe. 

The current measurements of the CMB by WMAP, the large catalogues of SDSS, the 
new cluster and the peculiar velocity catalogues motivate us to undertake considerations to 
explain features that are not expected in the concordance ACDM model. 

In the next chapter we derive the evolution equations of density-contrast for models with 
expansion, vorticity, acceleration and angular momentum within the corrected Ellis-Bruni 
formalism. A comparison is made with the standard formulas of isotropic and homogeneous 
spacetime. One can then evaluate integrated Sachs- Wolfe effect and RMS peculiar velocity 
field for any model. 

The last chapter is devoted for discussions, conclusions and suggestions for future work. 
Appendices contain detailed framework for Einstein-Cartan cosmology with field equations, 
propagation equations and a comparison of the corrected and standard Ellis-Bruni fiuid fiow 
approach for the density-contrast evolution. 

II. EVOLUTION EQUATIONS 

Our interest is the evolution of the Universe during the matter dominated epoch. The 
standard lore for the evolution of the density contrast, peculiar acceleration and peculiar 
velocity field gives the following equations [3, 8, 9]: 
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5{a, k) = 5{a)5t, 5{x) = / d^ke^'^ '^S 



The root mean square (RMS) values of the mass- density contrast and the pecuhar velocity 
field at a certain scale R and the redshift z using Gaussian window functions 3, S Q] are 
defined as: 



{6M/M)l^s{ci,R) = {{5M/Mf){a,R) = NV^^ j d'kW\k,R)5\a)\6^\\ (2) 
vlMsia.R) ^ {v'){a,R) = NV^' J d'kW^k, R)j^{ad^r\6^f, (3) 
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The aim is to study cosmological models with vorticity, acceleration and nonvanishing 
total angular momentum through torsion effects. We use the standard CDM power spectrum 



P{k) = defined in 

The growth function 6 (a) must be studied carefully. Our choice of covariant spacelike 
vectors within a fiuid fiow formalism differs from that of Ellis and Bruni jlo| : 



^ R\t)p-'h;V,p, (4) 
= R'{t)h;V,Q. 



ll| and their evolution equations result 



These vectors fulfil the Stewart- Walker lemma 
in a correct solution for a density-contrast formed from their scalar invariants. This is not 
the case for the standard Ellis-Bruni covariant vectors. 

The detailed derivation of the equations one can find in the Appendix A while a compar- 
ison between two fiuid-fiow approaches is in the Appendix B. 

The resulting second order coupled equations in our corrected scheme are given by (matter 
dominated regime): 
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V,. 



;^e5i + uy + Col + ^i)(^0^A - u^^a'-V, -V^- Vpia':^ + cu'',) + Gi^^a,) 

2 



^ A 



0, 



(5) 



We use the Zeljdovic model to describe the evolution of the total angular momentum of 



the Universe at smal 



galaxies and clusters 



redshifts assuming the surplus of right-handed over left-handed 



L{t) oc a^dp / d^q{q- f) x [V$o(g) - V$o(g)], 
L{Umverse) oc [n{right) — n{left)]L{t), 
Q = torsion oc L{Universe) =^ Q{a) = Q^a^^^"^ for z = a^^ — 1 < = 4, (6) 
where put arbitrarily, otherwise Q{z > z^r) = 0. 

The Einstein-Cartan equations remain unaltered with respect to the functional form of 
the time dependence of torsion (see Appendix A). 

We have to factorize density-contrast on the space and time dependent parts. One can 
achieve this goal transforming the evolution equations for covariant vectors to the local 
Lorentzian frame by tetrads: 

rf = diag{+l, -1, -1, -1), p,i^ = 0, 1, 2, 3, a,b = 6, 1, 2, 3. 
In the Appendix A one can find evaluated coefficients for the following spacetime metric: 



ds^ = dt" - R\t)[dx' + (1 - A2(t))e2'""rf|/2] - R^{t)dz^ - 2R{t)\{t)e^''dydt, (7) 

m = const. 
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It is easy to verify that in the Friedmann hmit one recovers the standard form of the 
density-contrast [l3| (see Appendix A): 



H(a) r . _3r ^0 13 



Ka)-^ldaa-^[j^r, (8) 
H{a) = Ho{nma-' + n^y/\ 

One can evaluate properly normalized peculiar velocities and integrated Sachs- Wolfe effect 



for various cosmological models ISj: 



-nm-^{5{a)/a), r = C daa-^H-\a), 5{a = 1) = 1. 



III. RESULTS AND DISCUSSION 

Supplied with all the necessary equations we can evaluate properly normalized peculiar 
velocities and integrated Sachs- Wolfe effect. We use the same normalization for all cosmo- 
logical models: 



{5M/M)RMs{.a = 1,R= lOh'^pc) = 1. 

Let us fix relevant parameters of the models in Table I (unit Hq = 1 is used for parameters 
m, Ao and Qq] EdS=Einstein-de Sitter, EC=Einstein-Cartan, ACZ)M=concordance model, 
= 1 — ^m, A = Ao-R~^, Q = QoR^^^'^, Rq = Hq^ ). 

The formation of small-scale structures and the age of the Universe can be explained 
with a larger mass-density and smaller Hubble constant (see Table II). This statement 
is valid if we assume that the total angular momentum of the Universe at low redshifts, 
acting through torsion terms, is much smaller than mass-density terms. The evolution 
with large torsion terms must contain feedback from matter to the background geometry, 
changing substantially its expansion and vorticity. It is possible to utilize this approach 
within N-body simulations. Thus, we limit numerical evaluations in this paper to small 
torsion contributions. 
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h 


m 


Ao 


Qo 


ACDM 


0.3 


0.7 











EdS 


1 


0.5 











ECl 


2 


0.4 


0.03 


0.067 


-0.2 


EC2 


2 


0.4 











ECS 


2 


0.4 


0.3 


0.67 






TABLE I: Model parameters 





ACDM 


EdS 


ECl 


EC2 


TuiGyi) 


13.77 


13.33 


13.14 


13.09 



TABLE II: Age of the Universe 



From Tables III- VI it is clear that only models with large mass-density can enhance 
large-scale peculiar velocities observed in the analyses with combined cluster and WMAP 

nn 

data p^, 115[ . Rather small amounts of vorticity, acceleration or torsion do not essentially 
influence RMS velocities. However, components of a density-contrast T^i^i = 1,2,3 are not 
equal. We assume the standard scaling of the vorticity A(a) = AqO^^ =^ lo oc (see 
Apendix A for definitions). 

The density-contrasts normalized at zero-redshift do not depend on the initial cosmic scale 
factor, but a difference between components does depend. One can estimate the resulting 
angle between the axis of vorticity (z-axis) and the anisotropic bulk velocity. The angle 
depends on the initial redshift and the magnitude of the vorticity {uj(t) = ^m\{t)R(t)~^): 



a{initial) = 10 ^, a{final) = 1, model = ECl, —j- = -mAo = 10 ^ 

Ho 2 

(P?+I??)V2 

^ l{n{flow),h{axis)) = arctg — = 55.3°, 

1^3 1 

a{initial) = 10^^, a{final) = 1, model = ECl, but m = 0.15 : l{n{flow),h{axis)) = 57.25°, 
a{initial) = 10~^, a(final) = l, model = ECl : l{n{f low), n{axis)) = 72.3° . 

Since the metric describes rotation around z-axis, the dominant components of 
over(under)densities, the peculiar accelerations and the velocities are placed in the plain 



z \ R(Mpc) 


50 


187.5 


325 


462.5 


600 





511.67 


184.73 


111.75 


79.95 


62.21 


0.25 


515.49 


186.11 


112.58 


80.55 


62.67 


1 


467.16 


168.66 


102.03 


73.00 


56.80 



TABLE III: Vrmsikm/s) for the ACDM model 



z \ R(Mpc) 


50 


187.5 


325 


462.5 


600 





766.27 


245.82 


145.55 


103.27 


80.00 


0.25 


685.38 


218.87 


130.18 


92.37 


71.56 


1 


541.84 


173.82 


102.92 


73.03 


56.57 



TABLE IV: frms (fcfn/s) for the EdS model 



z \ R(Mpc) 


50 


187.5 


325 


462.5 


600 





1090.94 


335.81 


197.54 


139.83 


108.19 


0.25 


901.75 


277.58 


163.28 


115.58 


89.42 


1 


678.04 


208.71 


122.78 


86.91 


67.24 



TABLE V: Vrms{km/s) for the ECl model 



z \ R(Mpc) 


50 


187.5 


325 


462.5 


600 





1091.69 


336.04 


197.68 


139.92 


108.26 


0.25 


902.28 


277.74 


163.38 


115.65 


89.48 


1 


678.45 


208.84 


122.85 


86.96 


67.28 



TABLE VI: Vrms{km/s) for the EC2 model 



perpendicular to the axis of rotation (vorticity). We estimate the angle between the mea- 
sured directions of the axis of vorticity and the large-scale flows 



n{flow) = {l = 287°, 6 = 8°), h{axis) = (/ = 260°, b = 60°) 

=^ l{n{flow),h{axis)) = 53°. 

The reader can compare and visualize density contrasts and their derivatives for two 
crucial models (ACDM and EC) in Figs. 1 and 2. 

The integrated Sachs- Wolfe effect is negative for large mass-density models (EC) with 
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Density-contrasts 




■] I I I I I I I I I 

0.5 1 1.5 2 2.5 3 3.5 4 

redshift 



FIG. 1: Density-contrasts for the ACDM and ECl models. 



negative cosmological constant, while it is positive for ACDM, as can be seen in Fig. 3. The 
negative contribution can explain the small large-scale power of the CMB TT spectrum jirj l. 

It seems that the introduction of rotational degrees of freedom (torsion, spin, vorticity, 
angular momentum) is inevitable in order to understand and fit all observational data. Two 
scenarios emerge as viable resolutions: small Hubble constant with small amount of the total 
angular momentum of the Universe at present or possibly larger Hubble constant if the total 
angular momentum appears much larger. Torsion terms (linear and quadratic) always give 
a negative contribution to the effective mass-density, as it can be seen from Einstein-Cartan 
field equations (see Appendix A). Considerations with large angular momentum (torsion) of 
the Universe must include large-scale N-body numerical simulations. 
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FIG. 2: Gradients of the density-contrasts ^ for the ACDM and ECl models. 



IV. APPENDIX A 



We provide here a complete set of coventions, identities and equations for Einstein-Cartan 
theory. Let us start with definitions: 



g^" = v^Xrr\ Vab = dtag{+l, -1, -1, -1), 
/i, z/ = 0, 1, 2, 3, a,b = 6, i, 2, 3, 



dA _ q pA _ n -pA I -pA p/3 _ pA p/3 



Field equations and Ricci identities look as 
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FIG. 3: Integrated Sachs- Wolfe x function for the ACDM and ECl models. 



Qa = v^Q^, = [ah] = l/2(a6 - 6a), {ah) = l/2(a6 + ha) 

Conformal or Weyl tensor CaXni^ is defined as: 



RaX/xi^ — ^iQaixRxu — QavRxix — QXuRau + QXuRaf^) 
— ^RigandXu — gai^gXn) + CaXtiv 
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The energy-momentum tensor of the Weyssenhoff fluid is derived by Obukhov and Ko- 
rotky 0: 



T^u = -{p- A)9t.u + U^[U,{P +p)+ 2M"V;35i 



(10) 



The Ehlers-decomposition of the velocity-gradient can be written as 



u^u^ = 1, hf,^ = - u^u^, = u^VuU^,, = V , 



(11) 



(12) 



The vorticity is uniquely deflned by a variational principle (see Eq. (3.11) of [19,]): 

^ij = 'fii'^aVj)u°'gi,^, ij = i,2,3. 

The above two formulas for vorticity agree, but a formula for vorticity in (5.5) of has 
a wrong sign, as well as definitions of vorticity in 201] and 2l|]. Event ually , this confusion 
caused some wrong terms in the derivation of the evolution equations in |22(] , as it is pointed 



out in 



211 ] and later in 23] 



The standard procedure leads to the evolution equations (Frenkel condition employed 
n^Q% = 0): 



3 



ElaP = Cf^a^ipW^U^ 



We use the following identity: 
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^'^V,(/) - h-;e^^jy = a J + {u^% + a^>: + l^h^f'^V^f, (13) 
and the continuity equation in matter dominated epoch: 

7/'^V^p + pVX = (14) 
to derive the evolution equation for density-contrast Eq.(5). 

We evaluate the coefficients of the coupled evolution equations for components in the 
local Lorentzian frame X>a with the metric of Eq.(7) (time component "Dg can be set to zero 
because of the relation u"-Va = 0): 



p. + hfy + aikV'^ + = 0, i, j, k^l,2, 3, (15) 



ail = — -A kA U Ko — AA 

3 3 3^ 3 3i? 

AA + XQ - -A^— A— ^ -3 + A^ , 

3 3i? ^ 3 i?2 12^ r' 3R^ 

Im ■ 1 m R 

= 2^^+ 4^^+ 2^^^' = 

bu = 2^, 6i2 = -2Q + A^, 6i3 = 0, 

Ji, ri 



Im ■ 1 it 1 m R 

^ "2:R^"4:R^"r:R:R' 

4-2 1 5^2 1 ll^x; 
= --A -3«A-3Q - g'^P-y^AA 

lv^ 5m^, 5,9,i?,n 5, A, 9 li?, 

023 = 0, 621 = 2Q - A^, 622 = 2^, 623 = 0, 



031 = 0, a32 = 0, 

1^2 1a 2^2 1 5i? • I-- 

033 = — A kA. U KO AA AA 

3 3 3^3^37? 3 

2m^^ 2/^2 lx2"i^ 1^/ o x2n 

+ 3i;^«-3i?^ - 6^V-3ii<-^ + ^ )' 

&31 = 0, 632 = 0, 633 = 2—, 
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(^1 = C?3 = 0, 

ri li it 

+ A(4g2 + «(2A - p) - 6-gA + A(5A + 2A— + 9A-)) 

Jrt K ri 

+ |(A(3 + TA^) + A(17A2 + 2kK - 2Q'' -Kp+ ^XQ) 

ri ri 

+ a|(3 + 5A2))], 
for 3 : aij = -aji, bij = -bji. 

The symmetric parts of Einstein-Cartan equations are given as: 



(11) : 2(1 - A )- + (1 - A )— - - A - 5-AA - AA = /^A + Q - — 
(22) : 2- + (1 - 3A^)^ - ^ - 2-AA = «A + 0^ - 

(33) : (1 - A')(2| + - (4 - A^)^ - S^AA - A^ - AA = «A + 
(6i) : ^(A- + -A) = gA + 2Ag + 3Ag-, 

(12) : i|(A + 2a|) = 0. 

In the hmit of small A^ << 1 we combine 66 and ii components to approximate the time 
gradients of the cosmic scale factor: 



^ = 3(«(p+A)-e^ + — + (-)^), 

R 1,2 . 1 4^, 5mXQ l.m,,, 

5= 2(3«A-3.p + -0^--^-3(-)=). (16) 



The age of the Universe follows immediately: 



Tu{Gyr) = — / — + Q^a ^ - -Q" + ^ + — 17 

A = Aqo"^, Q = Qofl"^''^; Ao, Qo evaluated in the unit Hq = 1. 
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Note that even a term linear in torsion Q is negative because mX > 0{mX < 0) implies 
Q < 0(Q > 0). 



The relation (3.31) of Ref. 19|, as the equation of motion for the angular momentum of 
the Zeljdovic model, is no more valid. 

Let us finally write the CDM spectrum used in numerical evaluations 



Ak 



(18) 



V. APPENDIX B 



In this appendix it is pointed out the difference between the standard Ellis-Bruni approach 
and its corrected version applied in this paper. Their variables are defined as 10|: 



= R{t)h:v,e. 



(19) 



The same procedure as in Appendix A results in the following equation for a density- 
contrast in the matter dominated epoch: 



+ {Vx + V,{a\ + u':x))QQ5l + u,a^ + u^^ + 



+ ^6^0^ + Ga^] = 0. 
The corresponding equations in the local Lorentzian frame are: 



(20) 



V- + hijV-- + aikV'^ + rf. = 0, i,j,k = 1, 2, 3, 



(21) 
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an 
fell 



R 

2 h — 

R^ R 



1 , ^ 



. Am 2 
2^^ ' 



Im ■ 7 i? 3 m R 

"ir^'Ir^^Trr' 



'r^ 



Til — 

-2g + A-, 6i3 



1 • m 7 -R 3 Rm 

= ~Tr^Ir^~Trr\ 

R ; o 1 ^ i R 

a22 = ^-\'-Q'- -f^P - 2AA- 



023 = 0, 621 = 2Q- A^, 622 = 4^, 623 = 0, 
031 = 0, 032 = 0, 033 = 2— + - 2^/^' 



^31 



0, &32 = 0, 633 = 4- 



R 

R' 



di = d^ = 0, 

d2 = R[2}? + \^(\+%) + {^f\{Q + l3X^) 

K K 

+ A(4Q2 + «:(2A - p) - 6-AQ 

ri 

+ A(5A + 2A^ + 9a|)) 

+ ^(A(3 + 7A2) + A(17A2 + 2kA - 2Q^ 
R 

- np + ^\Q + \^{3 + 5\'))l 
for i ^ j ■■ aij = -dji, hj = -bji. 

The authors in [lo| equahze components of their variables with a scalar density- 
contrast 6. This is possible in the Friedmann limes when all components are equal. However 
even then, the scalar quantity formed from their variables must be ad hoc multiplied by the 
cosmic scale factor to achieve the correct result: 



S oc R{t)[-V^V''Y/\ 
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Comparison of the two fluid-flow approaches 




-2.5 I ^ ^ ^ ^ 1 

2 4 6 8 10 

redshift 



FIG. 4: Comparison between the two fluid-flow approaches for the model ECS: 6 = [— P^T"^]^/^ 
and 6 = difference = {6- 6)/6, S{z = 0) = 6{z = 0) = 1. 

Our corrected variables V^, on the contrary, give immediately good and correct Fried- 
mann limes: 

d oc [-V^V]'/^ = [-VaV]'/^ (22) 
Let us stress that beyond Fricdmannian geometry two quantites are not equal: 



1/2 



hence we use throughout our paper corrected variables V^. In Fig. 4 the reader can find 
comparison between two formulas when the vorticity and the acceleration do not vanish. 
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